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Abstract We show that irradiation of the accretion disc 
by the white dwarf limits the occurrence of dwarf nova 
outbursts in post nova accretion discs. After the nova ex- 
plosion, the white dwarf has to cool for up to ^lOOyr 
- depending on the orbital period (i.e., disc size) and the 
temperature of the white dwarf after the nova-eruption - 
before the disc can begin producing dwarf nova out- 
bursts. During this time the inner disc is maintained in 
a hot, ionised state. Using these calculations, we inter- 
pret the long-term photometric variability of the post 
nova V446Her (Nova Her 1960) which shows today regu- 
lar dwarf nova outbursts. As the white dwarf in V446Her 
continues to cool over the next '^10-20yr, we predict an 
increase in the amplitude of outbursts and a decrease of 
the outburst frequency, because the decreasing irradiation 
of the accretion disc should allow an increasing annular 
extent of the accretion disc to participate in limit cycle 
oscillations. 

Key words: accretion, accretion discs - binaries: close - 
stars: individual: V446Her - novae, cataclysmic variables. 



1. Introduction 

Although nova eruptions and dwarf nova outbursts may 
resemble each other in some cases|^, the physical mecha- 
nisms behind the outbursts are entirely different: a nova 
eruption arises when hydrogen-rich material accreted onto 
the surface of a white dwarf ignites under degenerate con- 
ditions (see ^tarrfieldet al. 1998 for a review); dwarf novae 
outbursts are thought to result from thermal instabilities 
associated with hydrogen ionisation in an accretion disc 
(see [Cannizzo 1993 for a review and Ludwig et al. 1994 
for a detailed parameter study). 

It is clear that the class of cataclysmic variables (CVs), 
close binaries containing a white dwarf accreting from a 



^ In fact, a number of large amplitude/low outburst fre- 
quency dwarf novae are found in the lists of nova remnants, 
the most famous examp le being the short-period dwarf nova 
WZ Sge ([Puerbeck 1987|). 



less massive main-sequence companion, provides numer- 
ous potential novae progenitors (pre-novae). In fact, al- 
most all CVs should suffer nova eruptions repeatedly dur- 
ing their lives. To our knowledge, no system became a 
nova after it had been classified as a CV. However, a 
large number of nova remnants (post novae) have been 
found to be CVs after they attracted attention by erupt- 
ing (novae may arise, of course, on a white dwarf accret- 
ing hydrogen-rich material in other environments, e.g. in 
a symbiotic binary). Among the wide variety of known 
CV subtypes, most post novae were found to have rather 
high mass transfer rates M, and, thus, fall into the class of 
novalike variables^ Only in a small number of cases were 
quasiperiodic brightenings observed in post novae, identi- 
fying those systems to be dwarf novae ( Livio 198E ). 

In this paper, we analyse the circumstances under 
which a post nova will evolve into a dwarf nova follow- 
ing the nova eruption. In the next two sections we calcu- 
late the irradiation of the accretion disc by the hot white 
dwarf, and estimate the amount of time required for a 
post nova to cool to the point that it may begin exhibit- 
ing dwarf nova outbursts. We then discuss our results in 
the context of the post nova with the best observational 
coverage to test our predictions - V446 Her. 

2. White dwarf cooling in classical novae 

Even though a large part of the accreted hydrogen-rich 
material is ejected during the nova eruption, a significant 
amount of hydrogen is left over in the remaining envelope, 



which rapidly returns to hydrostatic equilibrium (Mac 
Donald 199^ ). The nova remnant experiences steady-state 



hydrogen shell burning until the nuclear fuel is exhausted, 
and is consequently heated to effective temperatures of 
several 10^ K. The turn-off times of the burning envelope 
which have been observed with ROSAT for a small number 
of novae are of the order of years, e.g. 1.8 yr in V1974 Cyg 



Note that the published mass transfer rates were deter- 
mined generally under the assumption that the disc luminos- 
ity arises from viscous dissipation only. However, irradiation by 
the hot white dwarf contributes significantly to the observed 
disc luminosity in young post novae. 
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Figure 1. Irradiation limits on the occurrence of dwarf 
nova outbursts as a function of time after nova outburst 
and orbital period. We use initial effective temperatures 
Twd = 5 X 10^ K (upper curves) and Twd = 3 x 10^ K 
(lower curves). The solid lines represent calculations in- 
cluding a boundary layer luminosity corresponding to 
-^wd = lO^^gs"^, whereas the dashed lines are calculated 
using only the luminosity of the cooling white dwarf. Note 
that t — corresponds to the end of the hydrogen shell 
burning phase, which might last for up to 10 yr after 
the actual nova eruption. 



( |Krautter et al. 1996| ) and 10 yr in GQMus ( [Shanley et al 
1995|). [Additional support for turn-off times on the order 
of years comes from monitoring of the ultraviolet lumi- 
nosity of a number of post no vae following their eruption 
([Gonzalez- Riestra ct al. 1998[) . 

Priahiili (198d) modelled the evolution of a classical 
nova through a complete cycle of accretion, outburst, mass 
loss, decline and resumed accretion. During the decline, 
the cooling of the white dwarf can be fit with a power-law 
of the form 



(1) 



where L denotes the luminosity of the white dwarf and 
t the time after the nova explosion. Prialnik's theoreti- 



cal model is confirmed by Somers & Naylor (1999), who is 



derive the cooling rate of the white dwarf using B band 
observations of the irradiated secondary in V1500Cygni 
and find 



Tuchman et al. (199C) and Mineshige et al. (1990) calcu- 
late the thermal equilibrium structure of externally irra- 
diated accretion disc annuli assuming that the irradiation 



flux is thermalized in the photosphere of the disc. King 



(1997) explains the UV-dclay in dwarf novae taking into 



account irradiation from the white dwarf which truncates 



the inner disc. Hameury et al. (1999) confirm this claim 
only for hot white dwarfs (T^d > 40 000 K) and find that 
the depletion of the inner disc must create several small 
outbursts between the main outbursts, contrary to ob- 
servations. Leach et al. (1999 ) explain the low states of 
VYScl stars as being due to the interplay of irradiation 
from hot white dwarfs and low mass transfer rates. 



Fo r com pleteness we note that the motivation of King 
(1997 ) and [Hameury ct al. (1999|) was to account for the 



■'problem of the UV delay" in dwarf novae, although the 
effects of irradiation of the inner disc are of general interest 
for th e operation of the accretion disc limit cycle mecha- 
nism. 3mak (1998[ ) critically re-examined the basis for the 
claim of a problem, and finds that there does not seem to 
be one. Previous studies which had claimed there to be a 
problem had utilised various simplifications in their time 
dependent models which turned out to be critical - such 
as the neglect of a variable outer boundary. The contrac- 
tion of the accretion disc during quiescence has the effect 
of increasing the local surface density in the outer disc, 
and thereby promoting disc instabilities which begin at 



large radii. These so-called "outside- in" (Cannizzo et al 



1986) or "type A" (Smak 1984) outbursts in which a nar 



row spike of enhanced surface density propagates from the 
initial site of the instability to the inner edge, are able to 
produce outbursts with the observed delay. 

In this work we calculate the structure of irradiated 
accretion discs in non-magnetic post novae, taking into 
account the time dependent flux of the white dwarf cooling 
down from the hydrogen shell burning phase, as well as the 
luminosity of the boundary layer resulting from accretion 
onto the white dwarf. The luminosity of the white dwarf 
is given by 



(t) = 47ri?wd^crr^d(i) oc r 



1.14 



(3) 



where i?wd and T^d denote the radius and the effective 
temperature of the white dwarf, respectively, and a is the 
Stefan-Boltzman constant. The boundary layer luminosity 



L 



BL 



G A'fwd M 



Rwd 



(4) 



where G is the gravitational constant and, as in Btehle & 
King (1999| ) and [Leach ct al. (1999| ) we take ^bl = 0.5 



3. Irradiation of the disc by the white dwarf 

External irradiation of the disc by the hot central source 
can suppress the thermal instabilities particularly in the 
inner disc regions. In the context of soft X-ray transients, 



Assuming that the boundary layer luminosity is radiated 
by the entire surface of the white dwarf and the disc is 
geometrically thin, the time-dependent flux Fi„{t) irradi- 
ating the disc at the radius R is given by 



-Firr(i) = (1-/3)- 



Lwd{t) 1 r 

2 [arcsm p ^ 

2na i?„d 



p(l-p2)^],(5) 
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(Adams et al. 1988; King 1997), where [3 is the albedo, p — enings that are reminiscent of dwarf nova outbursts (Livic 
i?wd/-fi and t = Q aX the end of the hydrogen shell burning 1989; Honeycutt et al. 1998a). The post novae for which 



phase. Throughout this paper we follow King (1997) in 
adopting /3 = 0.5. 

If the irradiation temperature (cT^'^j. = Fi„) exceeds 
Th ^ 6500 K the hydrogen in the disc is fully ionised in- 
dependent of the accretion rate ( ^an Paradijs 199*^ ) Thus, 
setting Tin- > Th at the outer edge of the disc (i? = i?out) 
gives a limit of the irradiation flux which suppresses the 
disc instability. Using standard equations ( Eggleton 1983| ; 
Frank et al. 1992) and assuming that the outer disc radius 
^out is 70 % of the primary's Roche lobe radius, we obtain 
the outer radius of the accretion disc as function of the or- 
bital Period P and the binary mass ratio q = M2/M„d- 
We used g = 0.5 for all calculations throughout the paper 
as we are mainly interested in systems above the period 
gap. 

We can now estimate how long disc instabilities in the 



this outburst behaviour is best documented are GKPer 
( = Nova Per 1901), which shows outbursts of ^ 3 mag 
every - 1000 d, and V446Her ( = Nova Her 1960), which 
shows outbursts with an amplitude of ^ 2.5 mag and a 
recurrence time of ^ 23d (Honeycutt et al. 1998b). 

The presence of outbursts in GKPer is not surpris- 
ing: The critical mass transfer rate below which outbursts 



are possible is Mcr — 10 ° g s i?: 



?2.6, 
10 



-0.87 



where i?io 



is the outer disc radius in units of 10 " cm and mi is 



accreti jn discs of post novae are suppressed due to irradi- 
ation, where we assume either Af = or M — lO^^gs^^, 



the primary mass in solar units (Cannizzo et al. 1988). 
The large disc size (Porb = 2d) for GKPer means that 
Mcr ~ 2 X 10^^ gs"\ which ensures that disc instabil- 
ities are possible even at very high mass transfer rates 
and almost independent of the mass transfer rate. Indeed, 
the outbursts observed in GK Per have been successfully 
modelled within the disc insta bility scenario ( Cannizzo fc 
Kenyon 198^ ; |Kim et al. 1992|) . 



as typi zal for dwarf novae above the period gap (for much 



higher accretion rates, the disc remains in a stable, hot 
state anyway). Fig. 1 shows the results for initial temper- 
atures of rwd(O) = 5 X lO^K and rwd(O) = 3 x 10^ K. 



By contrast, V446Her has Porb — 4.97 h (Thorstensen 
k Taylor 2600| ), which places it among a number of "nor 



mal" dwarf novae longward of the 2 — 3 hr period gap. The 
dwarf nova nature of V446 Her is confirmed by its spec- 



trum which is typical for this class ( Thorstensen & Taylor 



Apparently, disc instabilities, and, hence, dwarf nova out- 
bursts, should typically be suppressed for ~5-50yr, and 
up to ^100 yr if the white dwarf is strongly heated dur- 
ing the nova eruption. It is also clear that the size of the 
disc is a crucial parameter: the irradiation from the white 
dwarf can suppress the instability over the entire disc only 
for orbital periods Porb ^ 20 hr. The contribution of the 
boundary layer luminosity becomes important only after 
~ 5yr at the earliest. 

It is important to note that our results are lower lim- 
its on the time scale on which irradiation from the white 



2000 ) . The presence of dwarf nova outbursts in V446 Her 



dwarf nuppresses dwarf nova outbursts in post novae for (1 



indicates that the accretion rate cannot be too high. 

In the previous section we assumed that the disc in- 
stability is suppressed for Tin- > 6500 K, which is a rea- 
sonable but somewhat simplified statement. In order to 
get more detailed results, we re-derive the accretion rate 
above which the disc instability is suppressed, Mcr, from 
calculations of the vertical structure of irradiated accre- 
tion discs. Our computations are based on the vertical 
str ucture code written by one of us (JKC) and described 
i n [Cannizzo fc Wheeler (1984 ) and [Cannizzo fc Cameron 



the given disc albedo: (a) in order to produce significant 
dwarf nova outbursts, the disc instability has to affect on 
considerable parts of the outer disc and not only the outer 
edge (we focus on this subject in the next section) and (b) 
the disc may be flared and therefore intercept more flux 
from the white dwarf than described by Eq.(5), which as- 
sumes a flat disc. 



4. Dwarf novae among post novae? 

Accretion disc instabilities and their observational conse- 
quence, dwarf nova outbursts, can occur in a non-magnetic 
cataclysmic variable, if under the assumption of stationary 
accretion, the disc at its outer rim is colder than the ioni- 
sation temperature of hydrogen, Th. As described above, 
two factors determine the disc temperature at a given ra- 
dius: the mass transfer rate and the irradiation field from 



which we have modified to include the effects of 
external irradiation. Assuming that the external flux is 
thermalized in the surface layer, we set 



Cr = Pvisr + cr Tt,. 



(6) 



at the outer boundary, where Tg is the photospheric 
temperature, Pvisc the energy flux arising from viscous 
dissipation and Ti„ the irradiation temperature defined 
above (Hameury et al. 1999). We calculate the modi- 
fied vertical structure for A/^d = [0.5, 1.3] Mq, a = 
[0.005,. .., 1], P = [0.1, ...,10] X lO^^cm and T-„, = 
[0, 12] X 10'^ K. For the non-irradiated disc (i.e. Tirr = 
0), we obtain 



Mcr = 9.5 X 10l5gs-lp2^64^" 01 



-0.88 



(7) 



the white dwarf. 



This result is in good agreement with values derived previ- 
ously (Cannizzo et al. 1988 ; Ludwig et al. 1994 ; Hameury 



et al. 199^ ) 



While most post novae seem to have no outbursts, a 
small number of post novae show repetitive optical bright- 



Fig. 2 shows examples of the thermal equilibrium 
curves for different irradiation temperatures. Irradiation 
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has a significant influence on the ionisation state of the 
disc for log(Tiri.) > 3.5. The critical accretion rate, be- 
low which the disc becomes unstable, decreases with in- 
creasing irradiation until the external flux becomes strong 
enough to hold the disc in the hot stable state independent 
of the accretion rate. We find that the irradiation temper- 
ature TirrjS for which Met vanishes is a slightly decreasing 
function of i?io and a, but nearly independent of mi. 



7382 K a 



-0.07 



R 



-0.03 
10 



(8) 



The decrease of Mcr between the non-irradiated case 
(Eq. 7) and the suppression of the instability is approx- 
imately given by: 



5. Discussion 

Taken at face value, the decrease of the optical bright- 
ness observed in V446Her over several decades following 
the nova eruption and the onset of dwarf nova outbursts 
~ 1 — 3 decades after the nova eruption could be inter- 
preted simply as a decline of the accretion rate, being 
higher than the critical accretion rate for at least the first 
decade after the eruption, and being below the critical rate 
in the 1990's. 

Such a decline of the accretion rate in post novae is the 
hallmark of the "hibernation scenario" , which was invoked 



Mr: 



7.2 



rpj.2 

ii'r,s , 



(9) 



where Tirr,s and Mcr given above. Using the result of our 
parameter study, i.e. M^"(Mwd, -R, Tl^r, a), combined with 
Eq. (5) and assuming Twd(O) = 3 x 10^ K after the nova 
event, we find that the limit cycle instability can operate 
now (i.e. 40 yr after the 1960 nova) only between ^ 1.2 x 
10^° cm and Ront ^ 3.8 x 10^° cm in V446Her (Fig. 3). 
This might be enough to produce significant dwarf nova 
outbursts but the continuing cooling of the white dwarf 
should lead in the future to an increase of the outburst 
amplitude as well as to longer quiescence periods. 

Assuming that the mass transfer rate in V446Her is 
indeed low enough to permit a disc instability limit cycle, 
we cor mparc now our catimatoa of irradiation aupprcaaion 



by Shara et al. (1986, see Bhara 1989 for a review) to 
solve discrepancies between nova theory and observations 
and to account for the low observed space density of CVs 
when compared to the number of nova eruptions in our 
Galaxy or M31. In the hibernation scenario, irradiation 
of the secondary star by the hot white dwarf should keep 
the mass transfer rate high after the nova eruption for 
a limited period of time, after which the mass transfer 
rate decreases to very low values. The post novae should, 
hence, appear at first as a novalike variables (for many 
decades to a century) with stable hot accretion discs, and 
evolve thereafter into inconspicuous low M CVs. Thus, 
in the hibernation scenario it appears plausible that most 
"fresh" post novae have indeed high accretion rates, and 
that the oldest recovered novae are intrinsically very faint 
(e.g [Shara et al. 1985| ). 



On observational grounds, the hibernation scenario has 



of dwa rf nova outburata with the long term photometric 
history of the system. 

The pre-nova showed brightness fluctuations between 
V 18 — 15, which bear some resemblance to dwarf nova 
outbursts, even though the rise time of the single well sam- 
pled "flare" see ms too slow co m pared to a typ ical dwarf 
nova outburst (Robinson 1975). ^tienon (197l| ) reports 9 
photometric measurements of the post nova obtained be- 
tween 22 September 1968, and 26 September 1970, i.e. 
8-10 yr after the nova eruption, all of which show the sys- 
tem at TOpg 15.8. Finally, the systematic RoboScope 



been challenged, e.g. Naylor et al. (1992) and Weight et al 
(1994), but obtained also some support : |Vogt (199q ) a^ 



photometry of Honcycutt et al. (1998a), starting in 1990, 
shows V446Her at a quiescent magnitude of y ~ 18 with 
regular outbursts reaching V 15.5. It appears very un- 
likely that all of Stienon's measurements caught the sys- 
tem at the maximum of an outburst, and we believe that 
dwarf nova outbursts were suppressed in V446 Her for at 
least ten years after the nova eruption in 1960. Unfor- 
tunately, we do not know how long the hydrogen shell 
burning lasted in V446Her, and it is not clear at which 
point between 1970 and 1990 the dwarf nova outbursts re- 
sumed. Therefore, we can not give a quantitative estimate 

of rwd(o). 



Ducrbcck (1992| ) found a gradual decrease in the visual 
brightness during the final decline of novae and interpret 
this result as an evidence for a decrease of the mass trans- 
fer rate in these systems, consistent with the hibernation 
scenario. However, the latter point needs to be considered 
with some care: as already mentioned in the introduction, 
irradiation by the hot white dwarf significantly contributes 
to the observed disc brightness. As the white dwarf cools 
from the nova outburst, the decreasing amount of irradi- 
ation will result in a decrease of the observed disc bright- 
ness, and, hence, mimick a decrease of the accretion rate. 

A valuable alternative to the hibernation model that 
can explain at least the fact that only in a small number 
of cases dwarf nova outbursts were observed in post novae 
(and pre novae) is the mass transfer cycle proposed by 
pCing et al. (1995 ) in which the accretion-induced irradia- 
tion of the secondary drives a limit cycle with a period of 
~ 10^ — 10'' yr. In this cycle, a CV spends similar times in 
states of high and low accretion rates, appearing as a no- 
valike variable or as a dwarf nova. It is then clear that the 
probability that a CV turns into a nova is highest during 
the phase of high mass transfer, which naturally explains 
why most post novae are novalike variables. Of course, in 
some low M CVs the accreted envelope will also reach the 
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critical mass for a nova eruption, and V446 Her might be 
such a case. 

In this scenario, pre novae and post novae would have 
the same characteristics, which is in agreement with the 
analysis of Robinson (1975| ) , who found that in most novae 
the pre-eruption and post eruption magnitudes are very 
similar. 



6. Conclusion 

We have calculated the effect that irradiation by the 
heated white dwarf in a post nova has on the structure 
of an accretion disc. We find that, even if the accretion 
rate in such a system is low enough to permit disc insta- 
bilities, irradiation from the white dwarf suppresses dwarf 
nova outbursts for up to ~ 100 yr. In the case of V446 Her 
our calculations predict an increase of the outburst am- 
plitude and a decrease of the outburst frequency as the 
white dwarf keeps on cooling down. We encourage long- 
term monitoring of V446 Her in order to detect any evo- 
lution of the outburst pattern. 
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Figure 2. The vertical equilibrium curves for irradiation 
temperatures of Ti„ = [0, 2, 4, 6, 8, 10, 12] x 10^ K (right to 

left) and a white dwarf mass of M„d = 1-3 Mq. Assuming 
a relatively high a and at a large radius, the disc instability 
is suppressed for an irradiation temperature Tj^r = 8000 K 
(top panel), while it is not fully suprcsscd for a lower a 
(middle panel) or at a smaller radius (bottom panel) and 
the same irradiation temperature Tjrr = 8000 K. The ex- 
ternal irradiation temperature necessary to suppress the 
instability increases with decreasing radius and decreasing 
a (see Eq. (8)). Note the different scalings. 
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Figure 3. The critical accretion rate and the irradiation 
temperature as a function of time for different disc annuli 
{R = [3.5,3,2.5,2,1.5,1] x 10^° cm), a = 0.3 and mi = 
0.6. 



